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i-^ ■ Abstract 



^ ■ The effects of hght, long-hved gluinos on 2 ^ 2 processes at hadron colliders are 

examined. Such particles can mediate single squark resonant production via qg — > 
q ^ qg which would significantly modify the dijet data sample. We find that squark 
masses in the range 130 < rriq < 694, 595, 573 GeV are excluded for gluino masses 
of 0.4, 1.3, 5.0 GeV from existing UA2 and Tevatron data on dijet bump searches and 
angular distributions. Run II of the Tevatron has the capability of excluding this 
scenario for squark masses up to ~ 1 TeV. 
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Supersymmetry is a compelling candidate for physics beyond the Standard Model 
(SM) and has engrossed both the theoretical and experimental communities. Most of the 
attention has been focused on the minimal version of supersymmetry (MSSM), however, 
many other incarnations of supersymmetry could exist. In most cases these non- minimal 
models can significantly alter supersymmetric phenomenology and the associated search 
strategies, and hence all consequences of such models must be examined before regions of 
supersymmetric parameter space can be positively excluded. Here, we examine one such 
non- minimal case: the light gluino scenario. In some models it is natural ||l|] for gluinos to 
be much lighter than, e.g., squarks if they acquire their masses radiatively. While several 
experiments presently cast doubt on the existence of the low mass gluino window {rrig <, 5 
GeV), it has yet to be conclusively ruled out (or verified). In fact, the experimental bounds 
on this possibility are surprisingly spotty and controversial as evidenced by the continual 
debate in the literature |0. It is thus imperative to examine all implications of this hypothesis 
in order to quell this dispute. In this work, we investigate an additional data sample which 
provides strong constraints on the light gluino scenario, namely 2^2 processes at high 
energy hadron colliders. 

The window for a very light gluino was pointed out^] many years ago and its effects 
have since been analyzed in a variety of processes. A resurgence of interest in this scenario 
surfaced with the relatively recent observation Q that an apparent discrepancy between the 
value of as measured from jet production at SLD and LEP and that discerned from low energy 
data is resolved by the slower running of in the presence of light gluinos. However, recent 
compilations of various determinations of Og no longer show evidence of such a discrepancy, 
within the errors, but also claim that the precision of each individual measurement is such 
that any anomalous effect up to the ~ 5% level may not be perceived. The most noticeable 
consequence of this model is that the standard signals for gluino and squark production are 
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modified in the presence of liglit gluinos. Tlie bounds on the gluino mass, rrig > 144 — 224 
GeV from the Tevatron^] (with the range being due to the assumed relative sizes of the 
squark and gluino masses), are invalidated in this case as they depend on the fact that the 
g is short-lived and decays with the characteristic missing energy signature. Thus to be 
light, gluinos must be long-lived and appear to hadronize as jets. Since they are unable to 
appear as free particles, light gluinos will indeed form hadrons, with the bound states having 
longer lifetimes, and fragment in such a way as to mimic jets in a high energy detector]^. If 
kinematically allowed, the gluino hadrons will eventually decay into a final state containing 
jets -|-Xi5 where Xi is the lightest neutralino. The crucial ingredient for detection is then the 
ability of the final state Xi to pass the detector's missing energy cuts, which depends, amongst 
other things, on how the g hadron fragments. It has been estimated 0] that for rrig > 5 GeV 
the g would have been detected at UAl. However, as the gluino mass decreases, the missing 
energy signal disappears altogether. Standard squark searches are also nullified in this model 
as now the primary decay is g ^ qg, which again, escapes searches based on missing energy. 
In this case, the squark mass bounds are reduced to rrig > Mz/2, with the mass constraint 
being extended to 50 — 60 GeV from precision electroweak measurements at SLC/LEP0. 
We expect LEP II to strengthen the squark mass bound to ^ 80 — 85 GeV. 

We now discuss the results from a variety of light gluino searches. At present, the 
least controversial bound on light gluinos is from a search by CUSB|T^ for radiative T decays 
into bound states of gluinos. They exclude the mass range ~ 1.5 — 3.5 GeV (regardless of 
the gluino lifetime), where the lower limit is approximate due to questions! IT | concerning the 



validity of perturbative QCD in this regime. ARGUS [0 looked for secondary vertices from 
Xb ggg with subsequent decay of the gluino bound states and constrained a small region 
in the gluino mass - lifetime parameter space; these results, however, also suffer from 
perturbative QCD uncertainties as well as those from fragmentation effects. Beam dump 



experiments |]T3[ have looked for secondary vertices from the decay of g hadrons and appear 
to disfavor hght gluinos for restricted regions of the gluino hfetime, but these results depend 
on (i) assumptions on the production cross sections of the gluino hadrons, (ii) the value of 
the squark mass (iii) the interactions of the lightest color-singlet supersymmetric particle, 
Xi, with the detector, and (iv) g fragmentation effects and decay models. Searches for new 



neutral particles at Fermilab exclude |jT^ 2 < nig < 4 GeV for g lifetimes in excess of 10 s. 
Jet angular distributions of decays of the Z into four jets and precision measurements of the 
QCD structure constants Ca,f and Tp have been shown to be particularly sensitive to the 
existence of light gluinos [|l5|, but critically depend [|l^] on currently uncalculated higher order 



QCD corrections and hence no firm conclusions can presently be drawn. The Z boson can 
decay into 2 gluinos, however the branching fraction is small[|l^ {B ~ 0.06%), and would be 
hidden underneath ordinary QCD events. The detection of light gluinos at HERA, through 



their effect on deep inelastic structure functions ||18|| or via their production in the 3 + 1 jet 



photoproduction cross section |jT9|], have also been shown to be difficult. 



In this study, we examine the effects of light, long-lived gluinos on dijet production in 
hadronic collisions. One would expect the influence of light ^'s to be large in such processes 
since they contribute at leading order in perturbation theory. It has been shown [pO|] , however. 



that competing effects tend to suppress their impact on the single jet inclusive spectrum. 
Nonetheless, we find that the infiuence of resonant squark production from the subprocess 
qg q ^ qg should not be neglected as it greatly modifies the dijet mass spectrum and 
places strong constraints on the light gluino window. Our conclusions avoid some of the 
aforementioned difficulties in constraining this scenario, as non-perturbative QCD effects 
are negligible at the energies considered here and our results are insensitive to a long g 
lifetime. The essential ingredients of this model for our analysis are (i) the evolution of a.^ 
is modified by the inclusion of light gluinos in the QCD j3 function, (ii) long-lived gluinos in 
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the final state hadronize as jets, and (iii) light gluinos contribute a non-negligible partonic 
content of the proton. This introduces several new 2^2 parton scattering processes, as well 
as modifying the Altarelli-Parisi evolution of the parton densities. Global fits of structure 
functions which include a light gluino distribution have been performed ||18||, and it has been 



found that the NLO g parton distributions are roughly three (five) times larger than that of 
the strange quark at large (small) x for very light gluinos, rrig <^ 1.5 GeV, and carry ~ 5% 
of the proton's momentum fraction at large for rrig = 5 GeV. 

We now proceed with our calculation. All 2^2 subprocesses have been evaluated; 
they naturally fall into three categories, (i) those of the SM, qq —>■ qq , qq —>■ qq, gg , 
(19 ^ 19: and gg qq,gg, (ii) all SM initiated 2 — > 2 processes with final state gluinos, 
qq,99 99, and (iii) all gluino initiated processes, qg ^ qg , gg ^ gg, and gg gg,gg. 
Note that resonant squark production appears in the latter set. Higher order 2^3 processes, 
including the new reactions]^ which produce g+jet and thus yield 3 jet final states once the 
squark decays, have not been included. The mass of the light gluino has also been neglected 
in the evaluation of the subprocess cross sections as the results should not be sensitive to rrig 



at the energy scales considered here. The parton distributions ||18|| of Riickl and Vogt have 
been used for rrig < 1.5 GeV and those of Roberts and Stirling for rrig = 5 GeV. These values 
of the gluino mass avoid all of the experimental constraints detailed above. The change in the 
evolution of has been taken into account by fixing as{Mz) to the world average value p| 
and then running it to the relevant scale using the appropriate 2-loop (3 functions. We note 
that the 3-loop light g (3 functions have only recently been determined [P^] . 

In evaluating the squark resonance contribution to the cross section, we have used 
the narrow width approximation, which is valid for T/m <, 0.1 and hence is reliable in 
this case. We have included a 10% contribution to the squark width for potential non-dijet 
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decays, i.e., Vq = 1.1 x T{q — > qg). This is conservative as dijet decays will be by far the 
dominant mode. The 10% figure should cover the additional weak decays q —>■ qXi-4 aiid 
q — s> 5X^2 5 whichever are kinematically allowed, as they are expected to have small branching 
fractions of order < 1 — 2% each and hence are suppressed compared to the dijet mode. We 
note that monojet signals from squark production in this scenario have been previously 
analyzed p3|. We have also assumed that there are 5 degenerate squarks, with equal masses 



for the left- and right-handed states. Our results are not dependent on this assumption, 
however, as the contribution of each squark fiavor to the resonance peak is weighted by the 
corresponding quark's parton density. Hence this supposition does not simply result in an 
overall multiplicative factor to the cross section. In fact, the charm and bottom squarks have 
essentially negligible contributions to the resonance peak. 

Experimentally, the dijet system consists of the two jets with the highest transverse 
momentum in the event. In all cases, except where noted, we apply the cuts used by the 



CDF Collaboration!^ in their dijet analyses. This corresponds to pt^ > 20 GeV, jr^i 2I < 2, 
where 771^2 are the pseudorapidities of the two leading jets, and | cos 6**1 < 2/3, with 9* being 
the parton-parton scattering angle in the center of mass frame. Following CDF, we evaluate 
these processes at the scale /i = pt. In Fig. 1 we display the dijet invariant mass and single 
jet inclusive pt distributions for the cases of rriq = 300 and 500 GeV, taking rrig = 0.4, 1.3 and 
5 GeV corresponding to the dotted, dashed, and solid curves, respectively. In the case of the 
Pt distributions, we assume |?7i^2| < 0.5, /i = pr/2, and no angular cuts are applied. We see 
that the resonance peaks stand out for all values of the gluino mass. Note the degradation 
of the cross section as the g mass increases. 

We now evaluate the dijet resonance cross section and compare it to searches for 
dijet mass peaks from the single production of new particles performed by hadron collider 



experiments P4|, p5| , |26| . Figure 2 presents the single squark production cross section in the 
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nifj (UA2) 


rrig (CDF) 


(DO) 


0.4-1.3 
5.0 


130 to 195 - 220 
130 to 170 


220 to 475 
240 to 455 


310 to 590 
320 to 460 



Table 1: The squark mass regions in GeV excluded by the searches for dijet resonances by 
the UA2 (at 90% C.L.), CDF and DO (at 95% C.L.) Collaborations for an assumed gluino 



mass. 



dijet channel as a function of the squark mass for various values of nig. Also displayed in 
the figure (dotted curve) is the upper limit on the production of dijet 
UA2@ at 90% C.L., as well as both (b) CDF|g and (c) DO[0 at the 95% C.L. In the DO 
case the applied cuts are somewhat different than those employed by CDF: |?7i,2| < 1 and 
1^1 < 1-6. We see that the three experiments combine to exclude substantial regions of 
the light gluino parameter space. The ranges of the squark masses which are ruled out for 
each value of rrig are summarized in Table |1]. We do not expect the bounds to drastically 
improve as rrig ^ as the squark resonance cross section is not appreciably changing as the 
gluino mass decreases (once rrig < 1.5 GeV) as shown in Fig. 2. A short analysis shows that 
the cross section for massless gluinos is approximately 1.3(1.6) times larger than that for the 
case of rrig = 0.4 GeV at low(high) dijet invariant masses. 



Dijet angular distributions are a well known test of QCD and probe of new physics 



and have recently been measured at the Tevatron|2J, |26|. Ordinary QCD processes have 
large t- and w-channel poles and are thus peaked in the forward direction, whereas, resonant 
squark production in the light gluino model will have a fiat distribution due to the spin- 
nature of the squark. A convenient angular variable to use is x = exp(|?7i — 7/2!). For 
the case of 2 — 2 parton scattering, this is related to the center of mass scattering angle as 
X = (1 + 1 cos 6**1) /(I — I cos 6**1). x = 1 then corresponds to cos 6** = 90°. As is well-known 
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Figure 1: The (a) dijet invariant mass and (b) single jet inclusive px distributions for the 
2 — > 2 processes described in the text for the two cases rriq — 300, and 500 GeV. The gluino 
mass is taken to be 0.4, 1.3, and 5.0 GeV corresponding to the dotted, dashed, and solid 
curves, respectively. 
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Figure 2: Comparison of the single squark production cross section in the dijet channel as a 
function of the squark mass with rrig — 0.4, 0.7, 1.1, 1.3 and 5.0 GeV (straight curves, dotted, 
dashed, dash-dotted, square-dotted, and solid from top to bottom, respectively) with the 
upper bound for the production of new dijet mass resonances from (a) UA2 at 90% C.L., 
(b) CDF and (c) DO at 95% C.L. (dotted curves). 
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the advantage of the x variable is that it removes the apparent singularities associated with 
the t— and u— channel poles present in QCD. Thus da/dx shows greater sensitivity to new 
physics which does not possess such poles than does da /d cos 6*. To show the influence of the 
production of squark resonances on this distribution we display in Fig. 3 the ratio of da/dx 
calculated in the light gluino model to that of the SM, i.e., = {da / dx\g) / {da / dxlsAi) , for 



three dijet invariant mass bins (as chosen by CDF||2^) assuming a q resonance lies within 



each bin. In calculating the SM distributions, we employed the MRSA' parton densities! 28 
In all cases, we see that squark production leads to an enhancement in the distribution at 
low values of x compared to the SM. This would result in an increase in the dijet rate near 



90°. Comparison with the corresponding figures presented by CDF[23| shows that this rise 



in da/dx would be easily observable so that squarks with the masses chosen here could be 
excluded. 

We now make this procedure more rigorous in order to determine if the angular 
distributions can extend the excluded regions listed in Table |]. Following the procedure used 
by CDFp3, we employ the variable R{x) = N{x < 2.5)/N{2.5 < x < 5), which is the ratio 



of the number of dijet events in the two ranges of x, for the five mass bins 241 < Mjj < 300 , 
300 < Mjj < 400 , 400 < Mjj < 517, 517 < A% < 625, and 625 < Mjj GeV. This variable 
has the advantages that it is not very sensitive to variations in the parton densities, to the 
choice of renormalization scale {e.g., fi = pt versus Mjj), or to next-to-leading order QCD 
corrections, and that it characterizes the shape of the angular distribution in a mass bin with 
a single number. We have incorporated the systematic errors, as determined by CDF, as 
well as the statistical errors in our analysis. The systematic errors are highly correlated, and 



we have reconstructed the full covariance matrix according to the prescription in Ref. p4| 



We then calculate R{x) in each Mjj bin with rrig = 0.4, 1.3 and 5 GeV for squark masses 
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Figure 3: The ratio A^^ = {da / dx\g) / {da / dx\sM) as a function of x with gluino masses of 
0.4, 1.3, and 5 GcV, corresponding to the dotted, dashed, and sohd curves, respectively. The 
dijet invariant mass bins and assumed squark masses are as labeled. 
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in the range 160 — 800 GeV, and perform a fit to the CDF results using their data and 
correlation matrix. Following the usual analysis procedure, we find the minimum value 
of for a given value of rrig and then determine the excluded range of rriq by examining 
the distribution as a function of the squark mass. For definiteness we perform a LO 
calculation taking the scale fi = pt- Our results are presented in Fig. 4 for each assumed 
value of the gluino mass. Note that the x^ minima are generally found in the limit of very 
large squark masses. In all cases the distributions display a similar shape with 5 peaks 
which are associated with the 5 mass bins used by CDF and are due to the fact that the 
greatest sensitivity to a squark resonance occurs when it coincides in mass with the lower 
end of a given bin, i.e., when the squark cross section is maximum. To be more specific, 
when rUq is light (< 241 GeV) and outside the dijet mass region examined by CDF, the x^ 
is small but increases as the squark mass gets closer to the edge of lowest mass bin and then 
peaks once the bin is entered. The sensitivity then decreases as niq approaches the high end 
of the mass bin. As the value of rriq rises there is a general loss in sensitivity due to decrease 
in statistics and the corresponding increase in the size of the errors. 

This analysis excludes at the 95% C.L. the rriq ranges 151 — 694, 166 — 595, and 172 — 
573 GeV for 777,3=0.4, 1.3, and 5 GeV, respectively. It thus both extends and complements the 
constraints obtained from the dijet peak searches. Here, we might expect improvements on 
these constraints for rrig — > due to the increased enhancement in Ny^ at x = 0. Combining 
these results with the bounds from the resonance searches excludes squark masses in the 
range 130 < rriq < 694, 595, 573 GeV for gluino masses of 0.4, 1.3, 5.0 GeV. 

In summary, we have examined the constraints on models with light gluinos by using 
both the cross section and angular distribution for dijet events observed at hadron colliders. 
The critical observation is that a light gluino can act as a partonic component of the proton 
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Figure 4: ^ distributions as a function of the squark mass following the analysis described 
in the text, assuming gluino masses of (a) 0.4 GeV, (b) 1.3 GeV, and (c) 5 GeV. The dotted 
horizontal line represents the 95% C.L. bound in each case. 
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thus leading to the resonant production qg ^ q ^ qg, provided the q is sufficiently light. 
From our analysis, it would appear that the survival of the light gluino case requires either 
a light q in the ~ 70 — 130 GeV range, or a heavy q with rriq > 600 — 700 GeV. From studies 
of the physics capabilities at Run II of the Tevatron||29||, we anticipate that this future data 
will be able to exclude or verify this model for squark masses up to ~ 1 TeV. High energy 
hadron colliders may thus provide the best testing ground for this scenario. 
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